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Abstract 

This  paper  presents  the  concept  and  design  of  a two-layer  robust  control  system  for  suppression  of  combustion  instabilities 
over  a wide  range  of  operation.  The  control  law  synthesis  based  on  a multiple-time-scale  model  of  combustion  dynamics. 
Control  actions  on  the  fast-time  scale  are  provided  by  secondary  fuel  injection,  realized  as  modulation  of  the  primary  fuel  flow, 
and  are  gain-scheduled  according  to  the  variations  of  mean-flow  temperature  and  velocity  on  the  slow-time  scale.  A linear 
parameter  varying  (LPV)  Z,2  -gain  control  law  is  formulated  in  the  setting  of  differential  game  theory.  Simulation  experiments 

have  been  conducted  to  evaluate  the  control  law  under  wide-range  operation  of  a generic  combustor  in  tenns  of  the  trade-off 
among:  (i)  fuel  injection  rate  and  pressure  oscillation;  (ii)  transient  and  steady  responses;  and  (iii)  stability  robustness  and 
performance. 


1.  Introduction 

Unsteady  motions  in  confined  combustion  chambers  often 
result  from  mean-flow  dependent  interactions  between  heat 
release  and  acoustic  oscillations.  The  underlying  physical 
mechanisms  involve  a wide  range  of  time  scales  that  can  be 
broadly  classified  into  two  categories:  one  representing  mean- 
flow  motions  on  a slow  time-scale  and  the  other  representing 
acoustic  oscillations  on  a fast  time-scale.  The  disparity  of  time 
scales  allows  for  the  development  of  a hierarchically  structured 
control  law  for  regulating  the  various  processes  involved.  The 
overall  combustor  system  may  consist  of  two  interacting 
modules:  slow-time  supervisory  controller  and  fast-time  flame 
controller,  as  illustrated  in  Fig.l  [1],  The  slow-time  supervisory 
controller  in  the  outer  loop  attempts  to  optimize  the  overall 
combustor  performance  and  pollutant  emissions  based  on  the 
measurements  of  quasi-steady  variables,  and  to  identify  the  flame 
control  parameters  that  determine  the  operating  points  of  the  fast 
dynamic  phenomena  such  as  pressure  oscillations.  The  flame 
controller  in  the  inner  loop  is  responsible  for  suppression  of 
combustion  instability  and  for  prevention  of  lean  blowout  of  the 
flame.  When  coupled  with  a global  engine  control  system,  such  a 
two-layer  strategy  will  lead  to  an  energy-efficient  and  clean 
system  for  propulsion  and  power-generation  applications. 

As  part  of  their  research  on  feedback  control  of  combustion 
instabilities  in  various  propulsion  systems,  Yang  and  co-workers 
established  a series  of  control  designs  with  distributed  actuation 
[2-6].  The  recent  work  by  Hong  et  al.  [6]  was  based  on  the  Hx  - 
optimization,  which  guaranteed  robust  stability  and  performance 
within  specified  bounds  of  model  and  parameter  uncertainties  as 
well  as  sensor  noise  and  plant  disturbance.  The  nominal  system 
parameters,  however,  were  treated  as  time-invariant  constants, 
thereby  rendering  the  control  laws  valid  only  for  narrow  range  of 
operating  conditions.  The  present  paper  attempts  to  remedy  this 
deficiency  by  taking  into  account  the  temporal  variations  of 
mean-flow  dynamics  in  such  a unified  manner  that  the  resultant 
control  laws  function  effectively  over  a wide-range  of  operations. 
The  work  incorporates  an  -based  linear-parameter- varying 
(LPV)  control  [7,8]  for  modulating  combustion  dynamics.  Its 
major  contributions  beyond  the  previous  studies  on  narrow- range 
control  [6]  are 


• establishment  of  a two-time-scale  model  of  combustion 
dynamics  for  robust  control  systems  analysis  and  synthesis, 
and 

• development  of  an  observer-embedded  robust  LPV  control 
law  within  the  aforementioned  two- layer  architecture. 

The  paper  is  organized  in  six  sections  including  the 
introduction.  Section  2 develops  a two-time-scale  model  of 
combustion  dynamics  in  a liquid-fueled  propulsion  system.  The 
wide-range  robust  control  law  is  formulated  and  analyzed  in 
Sections  3 and  4,  respectively,  based  on  the  concept  of 
differential  game  theory.  Section  5 presents  the  results  of 
simulation  experiments  for  evaluating  the  system  performance 
under  different  control  laws.  The  paper  is  summarized  and 
concluded  in  Section  6 with  recommendations  for  future 
research. 

2.  A Two-Time-Scale  Model  of  Combustion  Dynamics 

The  disparate  time  scales  involved  in  the  various 
combustion  and  flow  processes  allows  for  the  decomposition  of 
the  conservation  laws  into  two  settings  governing  slow-time 
mean-flow  dynamics  and  fast-time  acoustic  dynamics, 
respectively.  The  fast-time  dynamics  is  parametrically  dependent 
on  the  slow-time  mean-flow  variables.  Using  a spatial-averaging 
procedure  equivalent  to  modal  analysis,  a family  of  linear  finite- 
dimensional  parameterized  models  is  derived  for  parameter- 
dependent  robust  feedback  control  design. 

The  equations  governing  the  dynamics  of  a liquid-fueled 
combustion  device  with  feedback  control  actuators  can  be 
succinctly  expressed  as  [2]: 


^■  + vg  -Vp  + pV- vg  =W  + I/VC 

(1) 

9v„ 

p-^-  + pvg  • Vv g + Vp  = F + FC 

(2) 

^+P'Vj+vg -Vp  = P+Pc 

(3) 

where  the  source  tenns  W , F , and  P represent  the  two-phase 
interactions  and  combustion  influence,  and  Wc  , Fc  and  Pc  the 
total  mass,  momentum  and  energy  control  inputs,  respectively. 

It  is  apparent  from  experimental  or  computational 
observation  that  motions  in  the  fluid  field  take  place  in  two  time- 
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scales:  one  representing  the  mean-flow  motion  with  a slower 
time-scale  and  larger  amplitude-order,  and  the  other  representing 
acoustic  motion  with  fast  time-scale  and  smaller  amplitude-order. 
Now,  the  purpose  is  to  decompose  each  single  conservation  law, 
described  in  Eqs.  (1),  (2)  and  (3),  into  two  time  scales  describing 
the  slow  behavior  of  mean-flow  motions  and  primary  combustion 
inputs,  and  fast  behavior  of  acoustic  oscillations  and  secondary- 
fuel  injection  into  the  combustion  chamber.  The  components  of 

the  fluid-variable  vector  x = [p  \g  p]T  in  the  following 

table  have  their  own  time-scales  and  amplitude-orders  where  the 
subscripts  ^ and  / represent  the  slow  and  fast  time  scales, 
respectively. 


Table  I Two  Time-scaled  Motions 


Time  scales 

Amplitude  Orders 

Mean-flow 
motion  x(ts) 

ts  =£st 

As  ~esAt 

Acoustic 

motion 

x'Of) 

*f  = eft 

Af  ~e/At 

In  the  table,  the  referenced  time-scale  and  amplitude-order  are 
denoted  as  t and  At  , respectively,  and  the  referenced  time  t 
here  is  taken  as  the  regular  time  used  in  the  physical  laws.  The 
time-scale  and  amplitude  order  are  described  by  the  pair  ( es,e  f ), 


with  £f  ~ ol(£, ) where 


p'lQ) 

e 


-»  0 as  6 — > 0 . 


Remark  1:  When  observing  events  in  the  fast-time  scale,  t * ~ t 
implying  that  £ r ~ 1 and  £s  » 1 . Similarly,  when  observing 


events  in  the  slow-time  scale,  ts  ~ t implying  that  e v ~ 1 and 


£ f « 1 . 

The  combustion  control  inputs  ( Wc  ,FC,PC)  are  also 

decomposed  as  two  components:  one  is  the  primary  component 
with  the  same  (i.e.,  slow)  time-scale  and  amplitude-order  as  those 
of  the  mean-flow  motions,  and  the  other  is  secondary  component 
with  the  same  (i.e.,  fast)  time-scale  and  amplitude-order  as  those 
of  acoustic  oscillations.  The  former  primarily  affects  the  mean- 
flow  field,  while  the  latter  is  a modulated  component  influencing 
the  acoustic  field.  Similar  treatments  are  applied  to  the  two- 
phase  interaction  variables  ( W,F,P ) of  the  uncontrolled 
combustion  process. 

The  universally  independent  variable  in  time-domain  is  the 
regular  time  t . However,  for  separation  of  the  local  and  global 
behavior,  the  fluid,  combustion,  and  two-phase  interaction 
variables  are  treated  in  the  form: 

x(t)  = x[tf(t),ts(t)}  (4) 

It  follows  from  the  above  equation  that  there  are  two  components 
in  x(tf,ts):  x'(tf)  is  only  a function  of  the  fast  time  tf, 

governed  by  the  acoustic  motion  and  the  other  x(ts)  is  only  a 


function  of  the  slow  time  ts  following  the  mean-flow  motion. 

The  two-time  scale  model  is  generated  as  follows: 

When  observing  an  event  in  two  different  time-scales,  we 


may  separate  the  time-derivative  terms 


d\„ 


dt  dt 


dp 


and 


dp_ 

dt 


Eqs.  (1),  (2)  and  (3)  in  corresponding  two  time-scales  in  the 
following  format,  say: 


avg[r,t/(?),ts(?)]_dvg  dts  dvgdtf_3vg  3v' 


dt 


dt..  dt  dtf  dt 


£s+——£f  (5) 

dtf 


Jl  f Ul  dts  Ul f 

Equation  (5)  implies  the  fluid  variables  p,  vg  and  p can  be 
linearly  expanded  as,  say: 

p(  r,  t)  = p(r,  ts)  + p '(r,  tf ) (6) 

As  the  slow-time  variables  are  related  to  system  performance,  the 
fast-time  variables  become  related  to  the  internal  stability.  If 
proper  slow-time  dynamics  is  pursued  in  operation  of  the  system, 
based  on  Eq.  (6)  and  the  magnitude  orders  indicated  in  the  above 
table,  those  tenns  vgVp  + pV-vg;  pyg-V\g+Vp; 

ypW • v g +vg  Vp  ; ( W,F,P );  and  ( WC,FC,PC ) in  Eqs.  (1),  (2) 

and  (3)  are  decomposed  into  two  equations  by  Taylor  series 
expansion,  say,  as: 

pV-  y g (r,  ts,tf)  = pV-vg  (r,  t s ) + p'(r,  tf  )V • vg  (r,  ts ) 

+ p(r,  ts  )V  • vg  (r,  tj ■ ) + higherorder  terns 
which  follows  the  following  general  form 


(7) 


f[x(tf,ts )]  = f[x(ts  )]  + - j-  v,/,  )-''V  / ) 


dx  I 

- higher  order  terms  • 


(8) 


Finally,  the  term  p- 


dvc 


dt 


in  Eq.  (2)  is  separated  based  on  Eqs.  (5) 


and  (6)  as: 

~ , d\  dv'„ 

= (p  + p)  (—£,+-—£ f) 


dt 


dts 


dt 


(9) 


/ 


which  is  approximated  to  the  first-order  accuracy  as: 


dVg 

dt 


~ dv 
:£.<P  — + £,P 
dt,. 


,dvc 


. 3v' 


ar+e'-p  at 


T1  0°) 


/ 


By  using  the  above  first-order  approximations  (5),  (7)  and 
(TO),  the  conservation  equations,  (T),  (2)  and  (3),  are  decomposed 
in  two  time  scales.  Collection  of  those  slow-time  variables  yields 

Slow-time  Conservation  Law: 


Mass 

J^-  + vg  (ts ) ■ Vp(ts ) + p(ts ) V • \g(t, ) 

= W(ts)  + Wc(ts) 

Momentum 

dvg(ts) 


£sp(ts) 


Energy 


dtL 

= F(ts)  + Fc(ts ) 
dp 


+ P(ts  )v  „ (ts ) ■ Vv  (ts ) + Vp(ts ) 


)V  ’ vg  (* s)  + Vg  Os  ) ■ Vp(f4  ) 
= P(ts)+Pc(ts) 

And  collections  of  fast-time  variables  yields 

Fast-time  Conservation  Law : 


Mass 

dp'(tf) 


-f  ■ 


dt 


f 


+ Vg  Os ) ■ Vp'O  f ) + y'gOf ) • Vp(ts ) 


+ pOs  )V  • v'g  O f ) + p'Of  )V  • Vg  Os ) 
= W'Of)+W'Of) 


(11) 


(12) 


(13) 


(14) 
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Momentum 

£/P(L)  '>y^0/  * + vP((f)  = ~P(ts)(\g(ts)-  Wg(tf) 

dt  j- 

+ Vg( tf)-V\g(ts))~  p'(tf)(es  1 (15) 

+ vg(t?)  ■ V\g(ts))  + F(tf)  + F'(tf) 


Energy 

e f + yp(ts  )V  • v'g  (tf)  = -yp\t  y)V  • \g(ts) 

Ot  f 

- yg(ts )Vp{tf)-Vg{tf)-Vp(ts)+P'(tf)+Pc{tf) 


The  set  of  conservation  equations  (11)  to  (13)  in  the  slow-time 
scale  governs  the  mean-flow  dynamics,  and  the  set  of 
conservation  equations  (14)  to  (16)  describes  the  mean-flow- 
dependent  acoustic  dynamics  in  the  fast-time  scale. 

Subtracting  the  fast-time  derivative  of  the  energy  equation 
(16)  from  the  spatial  derivative  of  the  momentum  equation  (15) 
yields: 


vV- 


a2{ts ) dzt / 


d2p'  , 

T~Tef  = 1 " 


hr.  = VF,  -- 


-2 


_1 dP' 

(ts  ) 


(17) 

(18) 


where  the  speed  of  sound  a is  a local  thermodynamic  property 
= JfRT  slowly  varying  in  an  ideal  gas 


_ dp 
given  by  a = — 

pP 

mixture.  The  function  hc  represents  the  high-frequency  control 


inputs,  and  h contains  the  driving  and  dissipation  terms  that 
depend  on  the  spatial  and  temporal  properties  of  mean-flow 
temperature  and  velocity. 

The  acoustic  dynamics  is  mainly  dependent  on  the  mean 
properties  vg  and  T that  determines  the  temporally  and 
spatially  compressible  properties  of  the  fluid  medium  in  the 
combustion  chamber.  Locally  and  temporally,  low  T implies 
reduced  speed  of  sound  and  consequently  the  dominant  natural 
frequencies  in  the  combustion  medium,  and  vice  versa.  The 
mean-flow  velocity  \ g represents  the  instantaneous  and  spatial 

speed  of  the  coordinate  frame  in  which  an  acoustic  mass  element 
travels,  providing  a driving/dissipated  force  due  to  an 

accelerated/decelerated  coordinate  frame.  Therefore,  T(r,ts) 
and  Vg  (r ,ts ) are  used  as  gain- scheduling  parameters  for  control 
of  mean-flow  dependent  acoustic  dynamics. 

Since  T(r,ts)  and  \g(r,ts)  are  continuously  varying  with 

respect  to  both  temporal  and  spatial  variables,  measurements  with 
finitely  many  point  sensors  may  cause  loss  of  relevant 
information.  Furthermore,  even  if  distributed  sensors  are  used  to 
measure  the  spatial  distribution  of  T(r,ts)  and  \g(r,ts),  the 

gain- scheduling  control  law  is  likely  to  be  excessively  time- 
consuming  and  hence  unsuitable  for  on-line  computation.  To 
synthesize  a control  law  that  is  executable  in  real  time  on  an 
inexpensive  platform,  the  following  approximation  is  made  for 
spatial  representation  of  the  gain- scheduling  variables  T(r,ts) 
and  Vg(r ,ts) : 


The  mean  gas  velocity  vg  is  first  separated  as  the  product 
of  the  temporal  variable  ts  and  the  spatial  (vector)  variable  r as: 


\g(ts,r)»\t(ts)<t>v(r) 


(19) 


T(ts,r)»Tr(ts)<l>T(r) 
where  the  scalar-valued  functions  0V( r)  and  <j>T( r)  are  non- 
negative at  each  spatial  point  r . The  mass  flow  rate  of  gas 
through  a combustor  of  cross  section  A( r)  is  given  by: 
m(t)  = A(r)\g(ts,r)p(ts,r) 

implying  that  the  mean  gas  density  p is  inversely  proportional  to 
^(r)Vg.  Here,  the  spatial  variable  r only  includes  the 

longitudinal  direction  after  spatial  averaging  over  the  cross- 
section. 

Since  the  mean  pressure  p has  been  assumed  to  be 
independent  of  the  spatial  variable  r in  high-pressure 
combustors,  the  ideal  gas  law  dictates  that  T is  inversely 
proportional  to  p . This,  in  turn,  implies  that  T is  directly 
proportional  to  A(r)  v „ as: 

f(tSfr)  = J^-A(r)vg(tJ,r)  (20) 

Rm(t) 

Using  Eq.  (19)  in  Eq.  (20)  yields: 

Tt(fs)  = ^-yt(ts) 

s Rm(t)  * (21) 

4>,  (v)  = A(v)<t>v  (r) 

For  a combustor  of  uniform  cross-section  (i.e.,  A(r)=A  for 
all  r ),  the  constant  term  A can  be  absorbed  in  the  temporal  part 
T,  (t)  so  that  Eq.  (2 1 ) is  modified  as: 


(22) 


Rm(ts) 

0/(r)  = 0v(r) 

Therefore,  for  uniform  cross-section  of  the  combustor,  Eq.  (19) 
can  be  rewritten  as: 

Vg((„r)  = vt(ts)<l>(r) 

T{ts,r)»Tt{ts)<p(r) 
where  the  function  0(r)  is  normalized  over  the  total  volume 
Vcomt  of  the  combustor,  i.e., 


(23) 


irJJI 


0(r  )dV  = 1 


(24) 


E comb 

In  the  sequel,  Eq.  (23)  will  be  used  for  modeling  combustion 
dynamics  and  control  systems  analysis  based  on  the  gain- 
scheduling variables  vt(ts)  and  Tt(ts)  that  are  available  from 
the  slow-time  scale  dynamics  of  the  combustion  control  system. 
The  scheduling  variables  vt  and  Tt  together  are  referred  to  as 
the  scheduling  vector: 

Tv, 

Tt 

The  speed  of  sound  a is  expressed  as 

a=at(t,)y[(j>(r)  (25) 

The  pressure  perturbation  can  be  modally  decomposed  as: 


p'(r,t  f)=^J)n(t  f)W  n(r) 


n= 1 


Space-time  separation  of  mean-flow  motions 


(26) 
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subject  to  specified  boundary  conditions  along  the  surface  of  the 
chamber.  In  this  coordinate  frame,  p has  a time-dependent 

component  rjk  along  the  direction  of  the  k,h  unit  vector  of  an 
orthonormal  basis  V»},!=i  2 > which  completely  spans  the  space 

of  pressure  functions.  That  is,  the  set  V«i,=i2—  satisfies  the 
following  conditions 


Onh„o™,liW: 


= 8„ 


for  all 


integers  m and  n ; 
Completeness:  (cp,  y/n  ) = JJJ 


<py/ndv  = 0 Vn  , 


then 


<p  = 0 . 

• Admissibility:  y/„  (r)  satisfies  specified  boundary  conditions 
Vn 

where  (•,•)  defines  the  inner  product. 

Modal  analysis  of  acoustic  pressure,  Eq.  (26),  is  performed 
as  follows. 


'VndV 


<n  = JJJ  ^JlmWmWndV  = JJJ  P 

m= 

=JJJ 


- 
dt 


a2h,.dV 


(27) 


= ^J]m  + Fn  + U„ 

m- 1 

where 

F.-fffrsr 

rrr  (28) 

l/»=-JJ  )a%dv 

In  view  of  Eq.  (27)  we  introduce  the  following  integral 
relationship: 

JJJ  a2  (r,  ts  )VymV„dV  = -0)]8mn  Vm,u  (29) 

to  approximate  the  mode  shapes  }„=i  1 ...  and  natural 
frequencies  {(9, , }„=1 , ...corresponding  to  specified  boundary 
conditions.  By  completeness  property  of  the  orthonormal  set 
V„}„= Eq.  (29)  is  equivalent  to  the  following  algebraic 
relationship: 


a 1 (r,  t, ) V V„  (r)  + ®„2  (t,  )W„  (r)  = 0 

(30) 

A set  of  ordinary  differential  equations  is  derived  as  an 

approximation  of  Eqs.  (27)  and  (28)  based  on  the 
separation  of  mean-flow  conditions  in  Eq.  (23): 

space-time 

+ o)2„  (Tt  (t,  = Fn  (t,  ,tf)+un  (t, , tf ) 

(31) 

where 

F"  “-«,2JJj0(r)/'^„(r)t/U 

(32) 

(33) 

and  a>n(ts)  = a,(t,)kn  (34) 

where  co„  is  the  slowly  varying  frequency  of  the  nth  mode  shape, 
and  k„  is  the  corresponding  wave  number.  Based  on  the  Eqs. 
(25)  and  (30),  approximated  mode  shapes  are  derived  as: 
0(r)VV„  + k2y/n  =0  (35) 

subject  to  the  boundary  conditions,  where  the  spatially  dependent 
weighting  function  0(r)  puts  a smaller  weight  on  the  mode 
shapes  of  pressure  perturbation  in  the  low-temperature  region  and 
a larger  weight  in  the  high-temperature  region. 

The  fast-time  scale  control  command  for  secondary  fuel 
injection  is  superimposed  on  the  slow-time  scale  primary  fuel 
flow.  Control  actions  arising  from  the  distributed  combustion  of 
secondary  fuel  are  modeled  by  an  ensemble  of  point  actuators; 
the  output  of  each  actuator  is  determined  by  its  position,  the  local 
burning  characteristics  of  the  fuel,  and  the  time  delay  from  the 
instant  of  fuel  injection.  A detailed  derivation  of  the  actuator 
command  U n in  Eq.  (33)  is  presented  below  [2]: 


U„ 


Cv 


2^¥n(r,)U;(tf;v) 


, [r / -r,  (v)] 

u =bi(v)~ 


dt 


(36) 


(37) 


where  A Hc  is  the  heat  of  combustion  per  unit  fuel  mass  and  the 
injection  rate  of  secondary  fuel  is  represented  by  min  with  spatial 
distribution  density  b: , the  time  delay  r,  associated  with  the  zth 
point  actuator. 

Finally,  the  open-loop  plant  model  of  the  mean-flow 
dependent  acoustic  dynamics  (in  the  fast-time  scale)  becomes 


d rin-+m2n(f(ts)yt]n 


dt 


(38) 


'^^Dni('r(ts),Vg(ts))-~~ Eni(T(ts),Vg(ts))T]i(tj) : 


1=1 

RAHC 

C„ 


^Jk(nts),yg(tsWn( rk) 


dmin[tf  -Tk(\g(ts))] 
dt, 


k J 

The  fonnulation  described  above  is  considered  as  a wide- 
range  model  of  combustion  dynamics  and  provides  a framework 
for  accommodating  both  the  global  and  local  behavior  in  the  two 
time-scale  control  systems.  For  both  control  and  estimation,  a 
state-space  realization  of  the  combustion  dynamics  represents  a 
family  of  v -parameterized  models.  The  models  that  include 
uncertainties  and  effects  of  mean-flow  dependence  are  presented 
below: 


xp  =(Ap(v)+Ap(v))xp  +G](v)v+G2  (v)d  ^ 

y = Cxp  +9  ' ’ 

where  v is  the  scheduling  vector  consisting  of  mean-flow 
temperature  and  velocity,  d is  the  plant  disturbance  and  9 is  the 

sensor  noise.  The  state  vector  is  chosen  as  xp  = (g  g)T  with 


g=rj , and  ?7  = [t?i  t?2  ■■■,Vn]T  ■ Unsteady  heat  release  O'  and 
acoustic  pressure  p'  provides  effort  and  flow  information  in 
mechanical  systems,  where  their  phase  difference  detennines  the 
dissipated  or  liberated  energy.  For  the  nominal  linear  model,  the 
system  matrix  is  given  as: 
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Ap  = 


0 / 

- £HT,)-E(v ) —D(v) 


with  Q = diagfcof , co2 , • • • , mjf ) 


(40) 


The  input  vector  v(t)  associated  with  a set  of  point  actuators  is 
related  to  the  mass  injection  rate  of  the  secondary  fuel,  min , as: 
b\  (v)w,„  (t—Tj  (v)—Stj  (v)) 
h (y)min(t-T  2 (v)  -Sr  2 (v)) 


v(t)  = 


(41) 


\bM  (v)m,„ (t-rM  (_)- St M (v)) J 
where  5r(-  is  the  modeling  error  associated  with  the  time  delay  of 

the  control  input.  Equation  (41)  represents  the  actuator  dynamics 
from  the  excitation  of  the  point  actuator(s)  to  the  actual  energy 
release  of  control  fuel.  Since  the  mean  flow  is  slowly  varying 
compared  to  the  motion  of  the  secondary-fuel  injector,  the  mean- 
flow-dependent  transfer  function  of  actuator  dynamics  can  be 
approximated  as: 

A-Sv)=^(T)e”Tl(^  bi (v)e*2Ws  •••  bM(y)e~T^)sJ  (42) 
An  additive  operator  Ap(v),  that  is  dependent  on  the 

mean-flow  conditions,  represents  the  effects  of  model  and 
parametric  uncertainties  of  the  plant.  It  can  be  treated  as  a 
disturbance  to  the  plant,  ws  =Ap(xp),  to  realize  the  energy 

amplification  from  input  to  output  through  this  operator.  The 
global  behavior  of  the  operator  A p is  characterized  by  the  L2  - 

gain  as  follows: 

lA^Lo^v)  (43) 

A physical  interpretation  of  Eq.  (43)  is  that  individual  dynamics 
of  the  uncertainty  operator  Ap(\)  yields  to  the  following  energy 

amplified  relationship  from  its  input  xp  to  its  output  ws  as: 

J |w?  | rfrcj  (5;t2(v)||.v/,||  dt  VTe[0,°°)  (44) 

for  zero  initial  conditions. 


3.  Formulation  and  Analysis  of  the  Wide-Range  Robust 
Control  Law 

This  section  focuses  on  robust  performance  issues  of  the 
inner-layer  within  the  two-layer  wide-range  combustion  control 
system  as  an  extension  of  the  narrow-range  control  concept 
presented  in  the  previous  publication  [9],  The  acoustic  dynamics 
at  the  inner  layer  is  parametrically  dependent  on  mean-flow 
dynamics  at  the  outer-layer.  For  robustness,  the  inner-layer 
controller  is  synthesized  such  that  stability  is  guaranteed  over  a 
wide  operating  range  of  persistently  varying  mean-flow  variables 
and  the  associated  modeling  uncertainties.  For  performance,  the 
controller  should  have  the  capability  for  disturbance  rejection  and 
meeting  the  requirements  of  acoustic  energy  and  control  energy 
under  both  steady-state  and  transient  conditions.  The  role  of  the 
controller  design  in  this  section  is  to  simultaneously  satisfy  these 
two  objectives  for  persistently  varying  mean-flow  conditions 
over  a wide  range. 

3.1  Two-layer  Control  Structure 

In  Section  2,  the  fluid-combustion  dynamics  have  been 
decomposed  into  acoustic  dynamics  in  the  fast-time  scale  and 
mean-flow  dynamics  in  the  slow-time  scale,  in  which  the 
parameters  of  acoustic  dynamics  are  dependent  on  mean-flow 


temperature  and  velocity.  In  the  integrated  control  strategy,  the 
fluid-combustion  dynamics  are  transformed  into  a combination  of 
slowly  varying  trajectory  of  quasi-steady  equilibrium  points, 
defined  by  mean-flow  dynamics,  and  the  associated  linear 
perturbations  are  defined  by  the  mean-flow-dependent  acoustic 
dynamics  in  the  fast-time  scale.  While  the  temporal  trajectory  of 
mean-flow  variables  is  tracked  for  global  performance  by  the 
slow-time  controller,  the  fast-time  controller  locally  suppresses 
the  acoustic  motions  with  its  slowly  varying  parameters  being 
gain-scheduled  by  the  mean-flow  temperature  and  velocity.  The 
control  system  of  fluid-combustion  dynamics  is  therefore 
designed  to  have  a two-layer  hierarchical  structure.  The 
aforementioned  mean-flow  dynamics  represents  the  open-loop 
plant  in  the  outer-layer  whereas  the  acoustic  dynamics, 
parametrically  dependent  on  the  mean-flow  variables,  represent 
the  open-loop  plant  in  the  inner-layer.  The  temporal  trajectory  of 
mean  flow  variables  is  tracked  by  manipulating  the  primary-fuel 
flow  in  the  slow-time  scale,  while,  in  the  inner-layer,  the  acoustic 
dynamics  is  feedback-controlled  by  secondary-fuel  flow,  that  is 
realized  as  fast  modulation  of  the  primary-fuel  flow. 

Figure  1 and  figure  2 show  a schematic  view  of  the 
hierarchical  structure  of  the  two-layer  control  strategy.  The  inner- 
layer  controller  continuously  receives  gain-scheduled  signals 
from  the  outer  layer  to  update  its  parameters  that  are  related  to 
the  temporal  equilibrium  points  defined  by  the  mean-flow 
temperature  and  velocity,  in  response  to  the  change  of 
perturbation  dynamics  at  quasi-steady  equilibrium  points.  A 
slow-time  scale  mean-flow  controller  at  the  outer  layer  is 
responsible  for  global-performance  of  the  combustion  chamber. 
Its  design  is  well  known  in  industry  and  is  not  reported  in  this 
paper.  The  inner-layer  controller  processes  the  fast-time-scale 
signals  of  pressure  oscillations  and  slow-time-scale  signals  of 
gain-scheduled  variables  to  modulate  the  primary-fuel  flow.  The 
modulated  signal  then  manipulates  the  actuator(s)  of  the 
combustion  chamber  to  locally  suppress  perturbations  along  the 
temporal  trajectory  of  the  mean-flow  variables  for  wide-range 
control  of  the  combustion  dynamics. 

3.2  Analysis  of  the  Robust  Control  Law 

To  design  the  robust  controller  at  the  inner  layer,  a family  of 
v -parameterized  generalized  plants  P(s;  v)  is  first  constructed 
following  the  block  diagram  in  Figure  3.  It  includes  the  v- 
parameterized  nominal  plant,  actuator  dynamics,  modeling 
uncertainties  and  performance  requirements  for  synthesis  of  a 
wide-range  robust  control  law.  The  feedback  controller  K(s\\) 
is  derived  as  a function  of  the  gain-scheduling  vector  v for 
controller  adaptation  under  slowly  varying  operating  conditions. 
Figure  3 shows  the  internal  structure  of  the  v -parameterized 
generalized  plants.  The  metric  of  L2  -gain  facilitates 
transformation  of  the  coupled  performance-stability  robustness 
analysis  problem  into  a stability  robustness  problem  (Zhou  et  al., 
1996).  To  incorporate  the  disturbance  rejection  capabilities  for 
robust  performance,  the  control  objective  is  related  to  the  metric 
of  L2  gain,  indicating  energy  amplification  from  the  generic 
disturbance  w to  the  objective  variable  z . The  generic 
disturbance  w consists  of  disturbances  induced  by  plant 
uncertainties  ws , disturbances  induced  by  v -parameterized 
modeling  errors  of  the  time  delay  function  wT  , weighted  plant 
disturbances  wd  , and  weighted  sensor  noise  we  . The  objective 
variable  z consists  of  stability  variables  associated  with  plant 
uncertainty  zs  and  time-delay  errors  zr , and  performance 
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variables  associated  with  the  acoustic  pressure  response  zp  and 
the  secondary-fuel  injection  zM  . 

The  subsystems  in  the  generalized  plant  family  are  the 
nominal  plant  dynamics  specified  in  the  v -parameterized  state- 
space  realization,  (Ap , 6j , G2 , C) , shaping  filters  associated 


with  the  plant  disturbances  Wd  , sensor  noise  Wg  , performance 
weighting  functions,  Wp  and  Wu  , associated  with  pressure 

response  and  secondary-fuel  injection,  respectively,  and  the 
stability  weighting  function  Wx  associated  with  time-delay 


errors.  Note  that  the  plant  uncertainty  Ap  and  phase  uncertainty 


At  induced  by  time-delay  errors  are  not  included  because 
modeling  uncertainties  are  represented  as  uncertainty-induced 
disturbances,  ws  and  wT  . Two  shaping  filters  Wd  and  We  are 
incorporated  to  characterize  the  main  frequency  components  of 
the  plant  disturbances  and  sensor  noise.  The  performance 
weighting  functions  W p and  Wu  are  specified  for  the  desired 

frequency  responses  of  acoustic  motions  and  secondary-fuel 
injection,  respectively,  providing  the  trade-off  between  (high- 
frequency)  transient  and  (low-frequency)  steady-state  responses. 
Those  subsystems  are  derived  as  analytic  functions  of  the  gain- 
scheduling variable  v to  serve  as  the  database  of  the  LPV 
controller  synthesis. 

Fast-time  modulation  of  the  primary  fuel  is  realized  as  a 
secondary-fuel  injection.  In  modeling  the  secondary-fuel 
distribution,  the  possible  parametric  error  of  time  delay  <5r(v) 
causes  v -parameterized  phase  uncertainty  of  the  closed-loop. 

The  v -parameterized  transfer  function  e~ST^'>s  is  represented 
as  one  member  of  the  following  set: 

^l+ArO(u;v)lfrO'(u;v)  V®  Vv  : ||Ar(/®;v)||  < l}  (45) 


where  At . (y®;v)  accounts  for  the  v -parameterized  phase 
uncertainty  and  acts  as  the  magnitude  scaling  for  each  component 
of  v . The  v -parameterized  stability  weighting  function 
WT  ( ja>;~ ) specifies  allowable  joint  phase-magnitude  margins  in 
the  closed-loop  system  and  serves  as  a metric  of  robustness. 
Equation  (45)  is  equivalent  to: 

_ J < 7)|  Vft)  w (46) 


where  St( v)  = maxl^.  ( v)|  Vv  . For  the  second  source  of 

modeling  uncertainty,  the  v -parameterized  plant  uncertainty 
A p(s;v)  acts  as  an  internal  feedback  to  the  nominal  dynamics. 


It  accounts  for  modeling  inaccuracies  of  intrinsic  coupling 
between  flow  dynamics  and  combustion  responses.  The  bound  of 

A p(s;v)  is  specified  by  its  L2  gain  as  ||a ^ (^ , v)|  <Sp(v). 


Based  on  the  small  gain  theorem  (Zhou  et  al.,  1996),  a sufficient 
condition  of  robust  stability  of  the  closed-loop  system  is: 


H Vretfioo) 


(47) 


for  the  zero-state  initial  conditions  and  zs  being  equal  to 
Sp(v)x p . 

The  inequality  in  Eq.  (47)  implies  that  the  robust  controller 
can  stabilize  all  perturbed  plants  within  the  v -parameterized 
uncertainty  bounds  characterized  by  Sp(v)  and  5t(v)  without 
specifying  any  performance  requirements  that  provide  trade-off 
between  the  control  energy  and  the  acoustic  energy;  transient 


response  and  steady-state  response;  and  stability  robustness  and 
performance.  The  performance  weights  within  the  system 
bandwidth  are  included  in  the  v -parameterized  generalized 
plants  for  regulation  of  the  secondary-fuel  injection  and  pressure 
dynamics.  Shaping  filters  are  used  to  penalize  the  dominant 
frequency  components  of  plant  disturbances  and  sensor  noise  for 
effective  rejection.  The  ability  to  suppress  flow  oscillations  is 
quantified  by  a positive  quadratic  energy-like  function: 


Orthonormality  of  the  acoustic  mode  shapes  simplifies  Eq.  (48) 
as: 


H =M2  (49) 


The  combustion  chamber  tends  to  become  free  of  pressure 
oscillations  if  H in  Eq.  (49)  approaches  zero. 

To  achieve  small  steady-state  amplitude  and  short  settling 
time  of  pressure  oscillations,  a new  performance  variable  z is 

introduced  as: 


zp(s)  = Wp(s;v)t](s) 


(50) 


where  the  Laplace  Transform  of  a time-dependent  quantity  is 
denoted  by  "A";  and  the  weighting  function  Wp(s,v)  is  defined 


W = 

p 


W 


P 1 


0 


0 w 


p2 


0 


0 w 


pN 


(51) 


A procedure  for  selecting  the  performance-weighting  matrix 
Wp(s,v)  of  pressure  oscillations  is  outlined  below: 


The  Bode  plot  of  Wpi(s),i  = 1,2, — , JV  resembles  that  of  an 

integration  operator  within  the  bandwidth  around  each  natural 
frequency  ®(-  (v)  over  the  range  of  the  gain-scheduling  variable 


v . The  shape  of  each  W pi  (s)  outside  the  respective  bandwidth 
is  assigned  to  be  flat.  The  amplitude  of  each  W pi  (s)  represents 

the  extent  of  overshoot  of  the  amplitude  of  its  mode  shape  for 
trade-off  between  transient  and  steady  state  responses  in  the 
sense  that  a high  overshoot  usually  leads  to  a short  settling  time. 
Thus,  the  resulting  performance  variable  z p yields  small  steady- 
state  amplitude  and  short  settling  time  of  pressure  fluctuations. 
Note  that  the  bandwidth  of  Wpi(s)  could  be  altered  by  mean- 

flow-dependent  modeling  uncertainties.  The  nominal  bandwidth 
corresponds  to  the  natural  frequency  of  acoustic  oscillations.  To 
determine  the  possible  bandwidth  range  for  all  perturbed  plants, 
natural  frequencies  within  the  uncertainty  bound  of  v need  to  be 
calculated. 

The  other  performance  requirement  involves  secondary-fuel 
injection.  The  frequency  response  of  mass  flow  rate  of  secondary 
fuel  has  limited  bandwidth  and  overshoot  amplitude  due  to  the 
inertial  effects  of  the  fuel  flow.  To  this  end,  a performance 
weight  Wu(s,v)  is  incorporated  into  the  generalized  plant 


model,  and  a new  performance  variable  zu  is  defined  as: 

zu(s)  = Wu(s;w)u(s)  (52) 

where  u(s)  and  zu(s)  are  the  Laplace  transforms  of  min(t ) and 


z„(t),  respectively.  The  v -parameterized  performance 


weighting  Wu  (s,  v)  is  chosen  via  the  following  procedure: 
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The  Bode  plot  of  W„(s)  resembles  that  of  a derivative  operator 
within  the  frequency  range  of  interest.  A small  zu  implies  that 
the  high-frequency  component  of  secondary  fuel  flow  rate  mm 
has  been  filtered  out.  At  each  v , the  shape  of  Wu  outside  the 

bandwidth  is  assigned  to  be  flat.  The  bandwidth  range  is  altered 
by  uncertainties  in  the  (mean-flow-dependent)  scheduling  vector 
v . 

The  nominal  perfonnance  is  specified  by  the  relationship 
between  exogenous  inputs,  plant  disturbance  d and  sensor  noise 
9 , and  performance  variables,  zp  and  zu  . A shaping  filer  Wd 

is  included  in  the  generalized  plant  such  that  the  resulting  robust 
controller  focuses  on  rejecting  the  dominated  frequency 
components  of  plant  disturbance.  Therefore,  the  weighted  plant 
disturbance  wd  is  used  for  perfonnance  specifications,  instead 

of  the  disturbance  d itself,  where 

™<i(s)  = Wd(s)d(s)  (53) 

Similar  conditions  are  introduced  in  the  sensor  noise  9 as: 


wg(s)  = Wg(s)d(s)  (54) 

where  wg  is  the  weighted  sensor  noise  and  Wg  is  shaping  the 
sensor  noise  9 . 

A linear  parameter  varying  (LPV)  robust  controller  is 
designed  such  that  the  plant  disturbance  d and  sensor  noise  9 
have  minimum  effects  on  the  acoustic  motions  and  control 
actions  from  the  energy  perspectives.  The  nominal  perfonnance 
is  specified  by: 

f ||9(v)zp(0||2  +||z„(0||2<* 

T (55) 

< | (|K(0|2+|K(0|2)^  vre[o,~) 

where  q is  a v -parameterized  positive  scalar,  representing  the 
weighting  factor  of  acoustic  motion  at  each  v . As  q is  set 
larger,  then  the  response  of  acoustic  motions  is  more  emphasized 
than  that  of  control  actions,  and  vice  versa.  Furthermore, 
increasing  q implies  a better  rejection  ability  of  exogenous 
inputs. 

Based  on  Eq.  (47)  and  Eq.  (55),  a sufficient  condition  for 
robust  performance  (i.e.,  a combination  of  robust  stability  and 
nominal  performance)  of  the  control  system  is: 

-r  2 rT  2 

J \\z\  dt  < J HI  dt  vr  6 [0,  °o)  Vw  6 L2  [0,  T] ; 


w’s 

with  z = 

qzp 

; w = 

wd 

zu 

wg 

. "T  _ 

_WT_ 

(56) 


If  Eq.  (56)  holds,  then  the  controller  internally  stabilizes  the 
closed-loop  for  all  perturbed  plants  with  desired  performance, 
subject  to  an  uncertainty  bound. 

Following  the  internal  structure  in  Fig.  3,  the  generalized 
plant  family  P(s;v)  can  be  expressed,  similar  to  that  narrow- 
range  control  in  the  previous  publication  [9]  by  state-space 
realization  as: 


x = A(v)x  + (v)w +B  2 (v)u 


z = Cj  (v)x+Di2(v)u 


)’  = C 2 (v)x  + Z>2i(v)w 


Note  that  Eq.  (56)  is  conservative  because  it  ignores  the  structure 
of  the  uncertainly  robust  perfonnance  operator  A(s;v) , which 
consists  of  model  uncertainties  and  perfonnance-induced 
uncertainties.  The  generalized  plant  family  P(s;v)  in  Eq.  (57) 
can  be  compensated  to  make  the  requirements  of  robust 
perfonnance  specified  in  Eq.  (56)  less  conservative.  Figure  4 
shows  that  the  generalized  plant  family  is  compensated  by  two 

v -parameterized  compensators,  D ' (s,  v)  and  D(s,v).  The 
resulting  compensated  generalized  plant  represents  the  structure 
A without  altering  the  uncertainty  size  and  perfonnance 
requirements.  This  is  achieved  by  choosing  the  D -compensator 
as: 

D( s,  v) A(s,  v)  = A(s,  v)D(  s,  v)  Vv  (58) 

This  compensation  approach,  known  as  the  D-K  iteration 
(Packard  and  Doyle,  1993),  can  be  iterated  many  times. 


4.  Synthesis  of  the  Wide-Range  Robust  Control  Law 


The  inner-layer  control  law  is  synthesized  using  the  concept 
of  Linear  Parameter  Varying  (LPV)-  L2  -gain  methodology  that  is 
presented  in  this  section.  The  objective  of  the  LPV-  L2  -gain 

control  law  is  to  reject  the  effects  of  chamber  disturbances  and 
mean-flow  dependent  model  uncertainties  while  optimizing  the 
specified  performance.  The  LPV-  L2  -gain  controller  has  a H„  - 
structure  and  is  derived  as  a gain-scheduling  process 
accommodating  continuous  variations  of  the  mean-flow 
conditions.  The  major  challenge  is  to  synthesize  a family  of  v - 
parameterized  control  laws  that  guarantee  robust  perfonnance 
under  the  continuously  varying  mean-flow  conditions  over  a wide 
operating  range. 

The  LPV  L2  -gain  controller  consists  of  two  main 
components.  The  first  component  is  a v -parameterized  observer 
that  dynamically  estimates  the  states  of  the  v -parameterized 
generalized  plants  described  by  Eq.  (57).  The  estimator  structure 
is  not  formulated  as  an  exact  plant  state  observer  due  to  the 
presence  of  a calibration  term  wca)  (v)  that  protects  the  control 
system  from  being  excessively  sensitive  to  exogenous  inputs  and 
modeling  uncertainties,  entering  the  observer  state  equation.  The 
second  component  is  a state-feedback  control  gain  matrix,  which 
determines  the  control  action  based  on  the  estimated  state  x as 
seen  in  the  controller  configuration  of  Figure  5. 

The  remaining  task  lies  in  the  determination  of  the  v - 
parameterized  observer  matrix  L(w) , controller  matrix  K (v) , 
and  the  calibrated  term  wca/(v)  such  that  the  sufficient 
condition  for  robust  perfonnance  in  Eq.  (56)  holds.  Usage  of 
parameter-dependent  Lyapunov  functions  in  the  controller 
synthesis  leads  to  differential  matrix  inequalities  as  reported  by 
Flong  [9].  The  results  are  derived  based  on  the  following 
assumptions, 

Assumption #1:  The  parameter-dependent  matrices  Dl2  and 
D2l  are  nonnalized  as: 


£>i2(v) 


0 

/ 


D2,(v)D2ri(v)  - / ; 


(59) 


Assumption  #2:  Du  is  reduced  to  zero  by  loop  shifting  [10]; 
Assumption  #3:  D22  is  set  to  zero  without  loss  of  generality 

[io]; 


(57) 
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Assumption  #4:  The  parameter  vector  \(l)  is  a continuous 


function  belonging  to  a compact  subset  of  9?1  and  the  parameter 
variation  rate  \(t)  is  bounded  as  |v;(t)|  < /},•  Vi e 

Assumption  #5:  The  system  matrices  of  the  generalized  plant 
family,  (A,B],B2,C1,C2,Dl2,D2l  ),  are  continuous  functions 


of  v ; 

Assumption  #6:  The  control  matrix  K(\)  and  the  observer 
matrix  L(\)  are  continuously  differentiable  functions  of  v ; 

The  main  results  are  summarized  below  as  Theorem  1 below 
and  the  detailed  proof  has  been  presented  by  Hong  [9]. 

Theorem  1:  All  candidate  LPV  L2  -gain  controllers  for  y- 
performance  and  P -variation  are  synthesized  by  solving  the 
following  two  linearly  coupled  matrix  inequalities  for  X(v)  and 

Z(y): 


X-^+AxX-'-^P^-B^  B[ 


3 V, 

(I^ffpcpr' 

4 


- 1 
0 


0 

-i! 


<0(60) 


where 

and 


Ax=A-B2D{2Cl 


z-lAz+iz^+^±Pi^^Q^m-^4f2  ^ 4 

tfz4  -I  0 

Pi  0 -fl 

where  Az  = A + - B{D2\C2 

1 r 

F,  = — BTx 


<0  (61) 


Y 


f2  = ~(b2  x + dJ2C\  ) 

The  LPV  L2  -gain  control  law  in  Figure  5 is  derived  in  terms  of 
the  control  matrix 

K(v)  = -(Bl ( v)X( v)  + D{2  (v)C,  (v))  , (62) 

the  estimation  matrix  is 

L(v)  = Z(v)C2  (v) , (63) 

and  the  calibration  tenn  is: 

wmax  (*>  v)  = (Bi  (v)  - Z(v)C2  (v)D2 i (v))F,  (v)x  (64) 

The  resulting  dynamics  of  the  LPV  robust  controller  becomes 
x = [A  + 5,F,  + B2F2  - ZC2  (C2  + D2  ,F,  )]x + ZC2  y 


u = ~(B2  X + d[2C\  )x 
where  x is  the  estimated  state. 

The  differential  matrix  inequalities  in  Theorem  1 are 
discretized  by  a finite-element  approach  as  a series  of  solvable 
linear  matrix  inequalities  [9],  The  commercially  available 
package  of  Matlab-LMI  toolbox  that  numerically  solves  these 
inequalities  via  a convex  optimization  procedure  has  been  used 
for  controller  synthesis  in  this  paper. 

Remark3 : Theorem  1 enhances  the  theory  of  linear  parameter 
varying  (LPV)  control  in  the  following  two  areas: 


• The  features  of  the  LPV  L2  -gain  control  law  in  Theorem  1 
are  similar  to  those  of  the  existing  LTI  control  (Doyle  et 
al..  1989). 

• Unlike  the  LPV  control  law  of  Wu  et  al.  (1996),  the 
proposed  LPV  control  is  not  explicitly  dependent  on  the 
temporal  rates  of  the  scheduling  variables.  This  feature  is 
attractive  from  the  perspectives  of  real-time  implementation 
because  measurements  of  time  derivatives  of  the  scheduling 
variables  may  not  usually  be  available. 

5.  Simulation  Results  and  Discussion 

This  section  presents  a series  of  parametric  studies  and 
transient  simulations  of  a generic  combustor  model  to  validate 
the  wide-range  robust  control  methodology.  The  objective  is  to 
investigate  the  relationships  among  system  performance,  stability 
robustness,  and  plant  operation  based  on  the  transient  response 
and  frequency  spectrum  of  pressure  oscillations  under  different 
operation  scenarios.  The  relationships  among  the  uncertainty 
bounds  of  system  dynamics,  the  response  of  flow  oscillations  and 
control  actions,  and  the  allowable  variation  rates  of  mean-flow 
parameters  are  investigated  and  quantified. 

Table  II  lists  the  parameters  of  the  generic  combustor  that 
represents  typical  scenarios  encountered  in  practical  combustion 
chambers.  Table  III  lists  the  dimensionless  quantities  that  are 
needed  to  investigate  the  four  dominant  modes  of  longitudinal 
acoustic  oscillations  under  Hx  Control,  Gain-scheduled  Control, 

and  the  LPV- L2 -gain  Control.  Table  IV  lists  the  parameters 
needed  for  modeling  and  simulation. 

Two  bandwidth-limited  performances  weighting  functions 
W p and  Wu  have  been  fitted  as  continuous  functions  of  the 

scheduling  variables:  mean-flow  temperature  T and  velocity 
v g at  each  of  the  selected  modal  frequencies  ©],■•■,  (tip  . The 

shaping  filters  Wd  and  Wq  are  taken  as  unity  as  the  chamber 
disturbances  and  sensor  noise  are  assumed  to  be  white.  The 
system  robustness  is  represented  by  the  bound  dpps)  of  plant 
modeling  uncertainties.  The  system  performance  is  represented 
by  the  weighting  factor  of  pressure  response  q(y) . The  plant 
operation  is  characterized  by  the  bounds  of  absolute  values  of  the 
normalized  rates,  \pT\  and  |pv| , respectively,  of  the  scheduling 

variables  T and  vg  . Increasing  |pr|  and  |pv|  implies  that  the 

scheduling  variables  can  be  varied  faster  with  guaranteed  robust 
performance.  On  the  other  hand,  increasing  the  performance 
parameters  q implies  a more  stringent  requirement  of  system 
performance  in  terms  of  the  transient  oscillations  of  each  acoustic 
mode  and  the  ability  to  reject  chamber  acoustic  and  thermal 
noise.  Extensive  simulation  experiments  have  been  conducted  to 
realize  the  trade-off  relationship  among  system  robustness, 
performance  and  operation  of  the  control  system. 

The  dimensionless  bounds  of  the  variation  rates  \pT\  and 

N jointly  define  the  allowable  varying  rate  of  plant  operation 

from  the  perspectives  of  performance.  Figure  6 shows  the  trade- 
off between  system  performance  and  plant  operation.  For  the 
plant  uncertainty  bound  set  at  Sp  = 0.03  , the  family  of  curves  in 
Figure  6 serve  as  a constraint  under  which  the  LPV  robust 
controller  guarantees  system  stability  and  performance  in  the 
operating  range  of  dimensionless  scheduling  variables 
pr  e [0. 6,1.4]  and  pv  e [0.5,1. 5] . The  triple  ( | p j- 1 ,|p,,|  ,q  ) 
defines  the  design  point  which  is  evaluated  in  the  modeling 
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process  prior  to  the  controller  design.  For  example,  with  the 
maximum  variation  of  the  mean-temperature  rate  \pT\  being  0.1 
and  the  maximum  variation  of  the  mean-velocity  variation  rate 
|p,,|  being  0.016,  the  weighting  factor  of  system  performance  q 

should  be  less  than  or  equal  to  0.09  to  guarantee  robust 
performance.  This  is  a conservative  design  and,  in  the  simulation 
stage,  much  larger  q may  be  allowed.  A larger  q leads  to  a 
larger  risk,  since  it  implies  a larger  probability  that  the 
combustion  system  may  enter  into  a (potentially)  unstable  region 
beyond  what  is  specified  by  the  plant  uncertainty  model.  The 
insufficient  information  of  uncertainty  dynamics,  except  its  size 
specified  by  the  induced  L2  -norm,  prevents  precise 
quantification  of  the  region  of  risk.  Figure  6 qualitatively  reveals 
the  sensitivity  of  each  design  point.  Robust  performance  of  the 
closed-loop  system  is  sensitive  to  the  mean-velocity  variation  rate 
pv  because  the  L2  -gain  control  tends  to  maximize  the  entropy 

of  the  closed-loop  system.  The  effects  of  mean-velocity  rate  pv 
on  the  dissipation  or  driving  force  of  acoustic  dynamics  represent 
the  varying  condition  of  the  irreversible  process  in  the  acoustic 
field.  Since  the  entropy-oriented  scheme  is  sensitive  to  the 
controlled  irreversible  process,  the  L2  -gain  control  system  of 

acoustic  dynamics  is  sensitive  to  mean-velocity  variations. 
Figure  6 also  suggests  that  the  robust  performance  of  the  closed- 
loop  is  relatively  insensitive  to  the  mean-temperature  variation 
rate  pT  because  the  mean  temperature  determines  the  natural 
frequencies  of  the  acoustic  vibration.  This  phenomenon  also 
leads  to  the  fact  that  the  bandwidth  of  the  open-loop  system  is 
almost  the  same  as  that  of  the  closed-loop  system.  In  other  words, 
the  feedback  control  does  not  change  the  natural  frequencies  of 
the  system.  Therefore,  model  reduction  originally  performed  in 
the  open- loop  plant  remains  valid  after  the  feedback  control  is 
formed.  Consequently,  for  large  natural  frequencies,  the  control 
actions  would  require  fast  actuators. 

Figure  7 shows  that  the  robust  performance  of  the  closed- 
loop  system  is  not  very  sensitive  to  mean  temperature  variation 
rate  with  |pv|  being  set  to  0.01.  The  results  indicate  that 

robustness  of  the  closed-loop  system  becomes  sensitive  to  mean- 
velocity  variation  rate  only  in  the  region  of  high  performance. 
Figure  7 can  be  used  to  realize  a trade-off  between  robustness 
dp  and  performance  q . In  the  low-performance  region,  the 
trade-off  between  robustness  and  performance  is  almost  linear, 
i.e.,  the  increment  of  8p  is  approximately  proportional  to  the 
decrement  of  q . 

Two  plates  in  Figure  8 show  the  temporal  trajectories  of 
mean  velocity  and  mean  temperature  generated  from  CFD 
calculation.  In  the  simulation  of  the  closed-loop  system  for  the 
wide  range  of  mean- flow  conditions,  the  plant  uncertainty  bound 
is  set  to  be  5p  = 0.03  , time  delay  uncertainty  bound  Sr  = 3/ co0. 


plant  disturbance  intensity 


jq-1  ^o£vPo 

RAHC 


and  sensor  noise 


intensity  10  1 ' . Simulation  experiments  are  conducted 

RAHC 

under  three  different  controllers: 

• Controller#!:  Single  control; 


• Controller  #2:  Gain-scheduled  H „ control; 


• Controller  #3:  LPV  L2  -gain  control. 


Figure  9 shows  the  system  response  using  a single 
controller  for  the  entire  range  of  pT  e [0.6, 1.4]  and 
pv  e [0.5, 1.4]  where  the  system  response  exhibits  instability  and 
hence  the  performance  is  obviously  unacceptable.  Note  that  H M 

control  may  be  suitable  for  narrow-range  operation  as  discussed 
in  the  previous  work  (Hong  et.  al.,  1999),  but  not  for  wide-range 
operation  here.  Figure  10  shows  the  system  performance  under 
the  gain-scheduled  controller,  which  is  synthesized  without 
considering  the  variations  in  the  mean-flow  variables  T and  \g  . 

The  results  show  that  the  control  system  is  unstable.  This 
establishes  the  necessity  for  LPV  robust  control. 

Figures  11,  12,  13  and  14  show  the  responses  of  four  modes 
of  acoustic  oscillations  in  the  combustion  chamber  under  the 
LPV  robust  control.  The  total  acoustic  oscillations  are 
significantly  improved  since  the  controller  exhibits  the  ability  to 
reject  exogenous  inputs  on  the  first  two  modes  without  degrading 
the  responses  of  high-frequency  modes.  This  suggests  that  model 
reduction  has  not  caused  any  significant  loss  of  robustness  under 
the  LPV  control  for  wide-range  operation.  Therefore,  usage  of 
LPV  robust  control  may  not  require  a meticulous  consideration  of 
the  high-frequency  excitation  provided  that  the  scheduling 
variables  are  appropriately  rate-bounded.  It  is  also  shown  that 
the  closed-loop  system  is  asymptotically  stable  in  the  absence  of 
exogenous  inputs. 

The  embedded  observer  located  in  the  LPV  robust  controller 
is  capable  of  estimating  the  internal  states  of  acoustic  dynamics 
under  a wide  range  of  continuously  varying  mean-flow  condition 
as  seen  in  Figures  11,  12,  13  and  14.  The  embedded  observer 
provides  an  estimation  of  acoustic  vibration  of  each  mode.  The 
transient  and  steady-state  errors,  due  to  modeling  uncertainties 
and  exogenous  inputs,  are  shown  to  be  small.  This  implies  that 
the  calibration  terms  in  the  embedded  observer  are  able  to  keep 
the  system  away  from  the  high-risk  region,  and  yet  they  are  small 
enough  not  to  significantly  change  the  estimated  states.  In  the 
absence  of  modeling  uncertainties,  the  performance  of  the 
observer  embedded  in  L2  -gain  controller  approaches  that  of  the 
minimum  variance  state  estimator. 

Figure  15  shows  the  secondary- fuel  injection  in  the  LPV 
robust  control  system.  The  fast-time  modulation  does  not  cross 
4 % of  primary  fuel  flow  rate,  while  the  amplitude  of  acoustic 
pressure  oscillations  exceeds  10%  of  the  mean  pressure.  Figure 
16  shows  the  frequency  spectrum  of  the  normalized  pressure 
perturbations  before  and  after  the  controller  is  brought  in.  The 
responses  exhibit  two  frequency-scaled  components 
corresponding  to  two  time-scaled  motions. 

6.  Summary  and  Conclusions 

To  synthesize  a robust  control  law  for  wide-range 
operations,  a two-time-scale  model  of  combustion  dynamics  has 
been  formulated  by  decomposition  of  the  conservation  laws  and 
thermodynamic  state  relations  based  on  the  slowly  varying  mean- 
flow  and  fast  perturbed  conditions  of  the  process  variables.  The 
fast-time  scale  acoustic-flame  dynamics  are  parametrically 
dependent  on  the  mean-flow  temperature  and  velocity.  Control 
actions  on  the  fast-time  scale  are  provided  by  secondary  fuel 
injection  that  is  realized  as  fast  modulation  of  the  primary  fuel 
flow.  The  structure  of  control  system  follows  a two-layered 
hierarchy.  In  the  outer-layer,  mean-flow  dynamics  of  the 
combustion  process  is  controlled  in  the  slow-time  scale  by 
manipulating  the  primary-fuel  flow,  while  the  fast-time  acoustic 
controller  in  the  inner-layer  is  parametrically  gain-scheduled  by 
the  mean-flow  temperature  and  velocity  from  the  outer-layer. 
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In  the  fast-time  scale  inner-layer,  a robust  control  law  is 
formulated  based  on  the  concept  of  linear  parameter  varying 
(LPV)  L2  -gain  in  the  setting  of  differential  game  theory.  The 

LPV  L2  -gain  control  law  has  an  embedded  observer  that 

dynamically  estimates  the  plant  states.  Time  derivatives  of  the 
gain-scheduling  variables  need  not  be  computed  or  estimated  to 
realize  the  interactions  between  inner-layer  and  outer-layer 
segments  of  the  control  system. 

For  both  narrow-range  and  wide-range  conditions  of  the 
mean-flow  variables,  the  robust  control  law  is  designed  in  terms 
of  robustness,  performance  and  parameter  variation  rate. 
Procedures  for  selection  of  mean-low-dependent  and  frequency- 
dependent  weighting  functions  have  been  proposed  in  this  paper. 
Physics  of  the  thenno-acoustic  phenomena  have  a significant 
bearing  upon  the  selection  of  these  performance  weights  and  very 
little  iterations  are  needed  to  arrive  at  the  final  design  of  the 
robust  control  law.  In  addition,  a numerical  tool  for  synthesis  of 
the  LPV  robust  control  law  is  formulated  based  on  the  concept  of 
finite-element  analysis.  Parametric  studies  and  simulation 
experiments  have  been  performed  to  evaluate  the  robust  control 
system  in  terms  of  trade-off  between:  (i)  fuel  injection  and 
pressure  oscillations;  (ii)  transient  response  and  steady  response; 
and  (iii)  stability  robustness  and  performance. 

The  results  generated  in  this  paper  potentially  lead  to  a 
design  methodology  for  wide-range  robust  control  of  pressure 
oscillations  in  a generic  combustor.  The  findings  of  this  paper 
could  be  applied  to  a multitude  of  other  combustion  processes. 
Although  the  simulated  results  clearly  indicate  efficacy  of  the 
robust  controller  design  methodology,  several  issues  concerning 
implementation  of  robust  control  laws  need  to  be  addressed 
before  applications  to  operating  combustors. 
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Table  II.  Parameters  of  the  Generic  Combustion  Chamber 


Time  average  sound  speed  a0 

795.8  m/s 

Time  average  mean  temperature  T0 

1415  K 

Time  average  mean  pressure  p0 

0.5  Mpa 

Combustor  length  L 

1.0  m 

Combustor  diameter  dc 

0.6  m 

Specific  heat  ratio  f 

1.2 

Gas  constant  R 

373  J/(kg.K) 

Constant  volume  specific  heat  Cv 

1500  J/(kg.K) 

(MMH)  heat  of  combustion  ( A Hc ) 

5xl06  J/kg 

Average  fundamental  frequency  <W0  = 7ta0  / L 

1 000  rad/s 

Average  mean  velocity  V0 

300  m/s 

Average  mixed  density  p 

0.9473  kg/m  3 

Table  III  Dimensionless  Quantities 


Dimensionless  time  t' 

co0t 

Dimensionless  frequency  ft/ 

O)/CQ0 

Dimensionless  injected  rate  of  control-fuel  u 

-RAHC 

— — M 

(l)aC\p0 

Dimensionless  damping  ratio  D’ni 

DJ(O0 

Dimensionless  frequency  shifting  E'nj 

EJcol 

Dimensionless  gain-scheduling  mean  velocity  pv 

v,/vo 

Dimensionless  gain-scheduling  mean  temperature  pT 

TIT 

Dimensionless  pv  variation  rate  pv 

dpjdt' 

Dimensionless  pT  variation  rate  pT 

dpT  / dt' 

Table  IV  Parameters  Used  in  Simulation  and  Modeling 


Maximum  time  delay  error 

St  = 3/ft)0 

Plant  uncertainty  bound 

dp  = 0.03 

Weighting  factor 

<7  = 0.1 

White  plant  disturbance  intensity 

1 0 1 ^oE'vPo 
RAH,. 

White  sensor  noise  intensity 

IQ-1  ®oCvPo 
RAH, 

Initial  conditions 

an  impulse  0.01  p0/ sec. 
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Figure  1 Two-layer  Control  System  for  Modulating  Combustion  Dynamics 


Figure  2 Two-Layer  Robust  Control  System  for  Combustion  Dynamics 


Figure  3 Mean-flow  Parameterized.  Generalized  Plant  of  Combustion  Dynamics 


Figure  5 Synthesis  of  LPV  Robust  Control  System 


Figure  4 Compensated  Generalized  Plant  of  Combustion  Dynamics 
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Figure  6 Trade-off  between  Performance  and  Plant  Operation 


Dimensionless  Time 

Figure  8 Time  History  of  Mean  Velocity  and  Temperature 


Dimensionless  Time 

Figure  10  System  Response  under  Gain-Scheduled  Control 
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Figure  7 Trade-off  between  Robusmess  and  Performance 


Dimensionless  Time 


Figure  9 System  Response  under  H ^ Control 
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Dimensionless  Time 

Figure  1 1 Time  History  of  First-Mode  of  Pressure  Oscillation 
under  LPV-  Z,  - Gain  Control 
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Figure  12  Time  History  of  Second-Mode  of  Pressure  Oscillation 
under  LPV-  Ij  - Gain  Control 
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Figure  13  Time  History  of  Third-Mode  of  Pressure  Oscillation 
under  LPV-  Lj,  - Gain  Control 
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Figure  14  Time  History  of  Fourth-Mode  of  Pressure  Oscillation 
under  LPV-  Z,  - Gain  Control 


Figure  1 6 Frequency  Spectrum  of  Pressure  with  and  without 
LPV  L. 2 -gain  Control 


